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Abstract 

We have used molecular dynamics simulations to construct a microcanonical caloric curve for 
a 1415-atom Ni icosahedron. Prior to melting the Ni cluster exhibits static solid-liquid phase 
coexistence. Initially a partial icosahedral structure coexists with a partially wetting melt. However 
at energies very close to the melting point the icosahedral structure is replaced by a truncated 
decahedral structure which is almost fully wet by the melt. This structure remains until the 
cluster fully melts. The transition appears to be driven by a preference for the melt to wet the 
decahedral structure. 
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It is well known that in small fee clusters, noncrystalline structures such as icosahedra and 



decahedra can 
of the interior 



become stable as the total surface energy becomes comparable to the energy 
l|. Similarly, in small coexisting solid-liquid clusters [2, [J the energy of the 
solid-liquid interface makes a substantial contribution to the total energy of the cluster, and 
can strongly influence cluster properties ODD . For example, to avoid the cost of forming 



the interface, small clusters tend to avoid static coexistence. This leads to an S-bend in the 
microcanonical caloric curve p|, and the corresponding negative heat capacities observed in 
small sodium clusters M- In fact, in sufficiently small clusters static solid-liquid coexistence 
may not occur at all [2| allowing only a dynamic coexistence between the solid and liquid 
phase However, in larger clusters, just above this threshold, where the cost of forming 
the interface is significant but not prohibitive, it may be that the solid component assumes 
new structures in order to form a more favorable interface with the liquid component. This 
would have interesting consequences for the interpretation of premelting features in cluster 
caloric curves Q|, and may also offer new ways of controlling cluster structures. 

Indeed, the delicate balance of internal and surface energies that produce icosahedra 
and decahedra is further upset in a coexisting cluster by the fact that the solid region is 
reduced in size. While the fully solid cluster may be in the size range where one structure 
is stable, the smaller solid cluster coexisting with the liquid may be in a size range where 
another structure is stable. It is possible that a rich variety of structures may be observed 
in coexisting solid-liquid clusters as new energetic balances are established, and it is not 
difficult to imagine a sequence of structural transitions as the liquid fraction of the solid- 
liquid cluster increases. There may be some experimental evidence for such a scenario: Koga 



et al |19| have recently reported the observation of an icosahedral to decahedral transition 
in free gold particles as they were annealed close to their melting temperatures. It was 
suggested that this was a thermally activated transition driven by cluster thermodynamics. 
However, it is quite possible that the gold particles were in a coexisting solid-liquid state at 
temperatures just below the melting point, and in this case the solid-liquid interface may 
well have played a key role in the structural transition. 

In this letter we report the transition of the solid region in a coexisting solid-liquid nickel 
cluster from a partial icosahedron to a decahedron. We have applied molecular dynamics 
simulations to construct a caloric curve for a 1415-atom cluster using an embedded atom 



method (EAM) potential [11]. Initially we attempted to identify the stable structure at 
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zero temperature by comparing the potential energies of closed-shell truncated octahedra 
(including variants such as the (TO) - and (TO) + structures 1^|), cuboctahedra, icosahe- 
dra, and Marks k| and Ino decahedra . In figure Q we show the relaxed energies for 
the icosahedra sequence and the Marks decahedra sequence relative to a fit to energies of 
the truncated octahedra sequence. In Ni we find that below 2869 atoms the icosahedral 
structures are lowest in energy, between 2869 and 12298 atoms the decahedral sequence is 
more stable and above 12298 atoms the fee truncated octahedra become favorable. This is 
similar to the structural sequence predicted in Ref |lj] using the same potential, although 
there it was estimated that fee structures would become stable at sizes of 17000 atoms. From 
figure ^ we note that the icosahedron should be stable at a cluster size of 1415-atoms but 
that the 1389-atom Marks decahedron and the 1289-atom truncated octahedron are only 
slightly less favorable energetically. Thus a 1415-atom nickel cluster would seem to be a 
good candidate for observing structural transitions during solid-liquid coexistence. 

A previous molecular dynamics study of solid-liquid coexistence in a 1289-atom nickel 
truncated octahedron [15] found partial wetting of the solid by the melt (with the solid 
exposing a solid lens consisting of (lll)-facets) and then full wetting at temperatures closer 
to the melting point (with the solid completely covered by the melt). This study also 
used the same EAM potential as our simulations. Other molecular dynamics simulations of 
coexistence in metal clusters have studied copper [l(|, gold 12 1 and lead j^. 

The caloric curves were constructed in the constant energy (microcanonical) ensemble 
using the following procedure: at each fixed total energy the cluster was equilibrated for 
150000 time steps (where At = 2 fs) and then the kinetic energy was averaged over a further 
150000 steps to obtain a temperature. An energy increment of 0.6 meV/atom was used to 
adjust the total energy between simulations by a uniform scaling of the kinetic energy. This 
corresponds to a heating rate of 1 meV/atom/nanosecond. To identify and characterize solid- 
liquid coexistence, we follow Cleveland et al [la], using the bimodality of the distribution 
of diffusion coefficients to distinguish solid and liquid atoms. We have previously used this 
method in Ref to characterize the coexisting solid-liquid states in Pb clusters. 

Figure El shows the resulting caloric curve for the 1415-atom icosahedron. At total energies 
below E = —3.82 eV/atom the cluster is fully solid. The onset of coexistence occurs at 
energies above just this. For energies between E = —3.82 eV/atom and E = —3.77 eV/atom 
the structure of the cluster is that of an incomplete icosahedron partially wetted by a region 
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of melt: figure HU compares the icosahedron structure prior to heating with snapshots of the 
coexisting cluster at E — —3.78 eV/atom with the liquid atoms shaded in dark (center), 
and then removed (right). In the coexisting structure it is possible to distinguish 6 nearly 
complete tetrahedra of the 20 that make up a full icosahedron. This structure is very 
similar to the icosahedral solid-liquid structures seen in molecular dynamics simulations of 
lead clusters |9(. We note that the region of exposed (lll)-facets of the icosahedron do not 
resemble a lens in the sense of Ref |l5|, where it used to refer to a patch of exposed solid 
surrounded by liquid in a nearly wetted cluster. 

At approximately E = —3.77 eV/atom a transition occurs. This is visible in the caloric 
curve by the sudden increase in temperature at this energy indicating that the cluster has 
lowered its potential energy. Note that the snapshot in figure |3] is taken just prior to this 
transition at E = —3.78 eV/atom. In figure 0] we show snapshots from immediately after 
the transition. The top left snapshot shows exposed (lll)-facets of the solid; the top right 
snapshot, looking at the opposite side of the cluster, shows no exposed solid. The region 
of exposed crystal facets is much more lens-like than in the case of the icosahedron with 
the melt appear to wet the solid more completely, although analysis of the distribution of 
diffusion coefficients suggests the liquid fraction has dropped slightly. In the bottom pictures, 



a common neighbor analysis 



17j (CNA) has been performed using the index classification 



from Ref to identify the structure of the solid region. The liquid atoms have been 
removed to show only atoms with fee or hep symmetries. Two angles are shown: looking 
down the fivefold axis of the solid (bottom left) and looking side on to this axis down one 
of the twin planes (bottom right). The structure closely resembles that of a five-shell Marks 



decahedron 



131 ] . We conclude the that transition seen at E = —3.77 eV/atom in the caloric 



curve (figure |2J) is a transition of the solid from an incomplete icosahedron to a complete 
decahedron. 

Figure El shows a second caloric curve constructed by first heating from E = —3.83 
eV/atom to E — —3.76 eV/atom and then cooling back to E = —3.83 eV/atom at the same 
rate. The transition from icosahedron to decahedron occurs at E = —3.766 eV/atom, which 
is close to the transition in the first caloric curve (figure suggesting that this is a good 
estimate of the transition energy. However we note in figure |S] that upon cooling there is 
no transition back to the icosahedron. Presumably, the decahedron is kinetically trapped, a 
common occurrence in both simulations and real clusters The the decahedron is fully 
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solid at approximately E = —3.79 eV/atom. Note that the temperature of the decahedron 
falls below that of the icosahedron at E — —3.80 eV/atom; at this energy both clusters are 
solid indicating that the solid icosahedral structure has a lower potential energy than the 
decahedron consistent with the zero temperature calculations (figure P). 

Also shown in figure are the number of atoms in bulk fee positions (obtained via 
CNA analysis of the cluster structure) as the cluster is heated and cooled. During the 
heating phase the number of bulk fee atoms steadily declines (as the solid fraction of the 
icosahedron decreases) until the transition at E — —3.766 eV/atom. At this point the 
number of bulk fee atoms (Nf cc ) jumps sharply as the transition to the decahedral structure 
occurs. Upon cooling Nf cc increases (as the solid fraction increases) until coexistence ceases 
at about E = —3.79 eV/atom where it can be seen that Nf cc becomes relatively static. In 
a perfect 1389-atom Marks decahedron (that is a relaxed structure at zero temperature) 
CNA analysis counts Nf cc = 650 and in a perfect 1415-atom icosahedron CNA analysis 
counts Nf cc = 400. Both the decahedral and icosahedral structure approach these values by 
E = -3.83 eV/atom. 

To test whether the location of the transition depends on cooling rate we conducted longer 
4 ns (2 million steps) constant energy simulations of the icosahedron at several energies near 
the transition point in the caloric curve. In figure we show the time evolution of the 
temperature and the number of fee atoms in the E = —3.77 eV/atom simulation. Here 
the transition is seen to occur at t — 3.65 ns where an increase in Nf cc and a jump in 
the temperature is visible. We see no such transition in 4 ns simulations at E — —3.78 
eV/atom and E = —3.775 eV/atom. This again suggests that we have a good estimate of 
the transition energy. 

There are several known instances where solid clusters undergo structural transitions 
prior to melting. Small gold clusters have been seen in simulations to undergo transitions 
from truncated-octahedral structures, which are globally stable at zero temperature, to 
icosahedral structures at energies just below the melting point [20]. Similarly a study of 
icosahedral Morse clusters found complex surface reconstructions prior to melting [21j . In 
these cases, the transitions are driven by the thermodynamics of the solid phase rather than 
contact with a melt. However it is unlikely that similar thermodynamics are driving the 
transition seen here as the transition is accompanied by an decrease in potential energy of 
about 10 meV/atom (corresponding to the rise in temperature that can be seen in figure EJ. 
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This decrease in potential energy of the overall cluster comes from the contribution of surface 
atoms which drop on average by 40 meV/atom after the transition; it seems likely that this 
comes from an improvement in the solid-liquid interfacial energy associated with the wetting 
of the decahedral structure. We note that the interior atoms actually experience a net 
increase in potential energy of approximately 10 meV/atom. As remarked earlier, a previous 
study of phase coexistence in Ni clusters found a preference for fee cuboctahedral cluster to 
expose (lll)-facets (lenses) to the vacuum or vapour [3], suggesting a preference for the 
melt to wet the (lOO)-facets. Molten lead is known to have a preference to wet crystalline 
lead (lOO)-facets over the (lll)-facets 22!] and the structure of the coexisting icosahedron in 
: igure |3] strongly resembles the coexisting icosahedron structure seen in simulations of lead 
9|. It is possible that the higher-energy (lOO)-facets on the decahedron offer an improved 
solid-liquid interfacial energy over the (lll)-facets; this is hard to verify directly due to the 
dificulty in defining the solid-liquid interface in such a small system. 

At cluster sizes near 1400 atoms, the EAM potential [111] predicts that several structures 
possess very similar energies, including the icosahedron and the decahedron (with a size 
of 1389 atoms). Indeed the transition occurs between an incomplete icosahedron and a 

n 

decahedron. In simulations of Ni clusters at other sizes [23], we have observed static solid- 
liquid coexistence in a 923-atom icosahedron (but not in cluster sizes below this) and we did 
not observe this transition to a decahedral structure. 

We conclude that we have found a structural transition that occurs in the solid part of 
a coexisting solid-liquid nickel cluster modelled using an EAM potential. We believe this 
is the first suggestion that structural changes can occur in solid clusters to accommodate a 
coexisting melt. This effect may explain the recent observation of a icosahedral to decahedral 
transition seen in free gold particles 3] as they were annealed close to the melting point and 
it could provide an important kinetic mechanism for controlling the structure of nanoscale 
metal nanoparticles. 
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FIG. 1: Comparison of the relaxed zero temperature energies of icosahedra, Marks decahedra and 
truncated octahedra for nickel. Energies are given relative to a fit to the energies of the truncated 
octahedra sequence. 
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FIG. 2: Caloric curve for the 1415-atom icosahedron. The onset of coexistence occurs at an energy 
of approximately E = —3.81 eV/atom. There is a second transition at E = —3.77 eV/atom, 
followed finally by melting at E = —3.755 eV/atom. 




FIG. 3: The figures show icosahedron prior to heating (left) and the coexisting solid-liquid icosahe- 
dron at E = —3.78 eV/atom just prior to the transition (center and right). In the center the liquid 
atoms are shown in a darker shade. Note that the liquid only partially wets the solid, exposing 
(lll)-facets, whereas if the liquid completely wet the solid then there would be no solid exposed. 
On the right, the liquid atoms have been removed to show only the solid. 



FIG. 4: The plots show the coexisting solid-liquid decahedron at E = —3.77 eV/atom just after 
the transition. In the top two snapshots, taken from opposite viewpoints, the liquid atoms are 
shown in a darker shade. In the bottom snapshots, the liquid atoms have been removed and CNA 
analysis has been used to highlight the twin planes of the decahedron. On the bottom left the 
viewpoint is down the fivefold axis of the decahedron. On the right the picture is looking side on 
to one of the twin planes. 
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FIG. 5: The plot shows a second caloric curve where the cluster is heated (solid lines) until the 
icosahedral to decahedral transition and then cooled (dashed lines). Also shown is the correspond- 
ing the number of fee atoms (as calculated by CNA analysis) in the cluster. 
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FIG. 6: The plot shows the time evolution of the temperature and number of fee atoms in the 
cluster as it goes through the transition from icosahedron to decahedron at E = —3.77 (eV/atom). 



13 



